Green’s function for T-stress of semi-infinite plane crack

Green’s function for the T-stress near a crack tip is addressed with an analytic function method for a semi-infinite crack lying in an elastical, isotropic, and infinite plate. The cracked plate is loaded by a single inclined concentrated force at an interior point. The complex potentials are obtained based on a superposition principle, which provide the solutions to the plane problems of elasticity. The regular parts of the potentials are extracted in an asymptotic analysis. Based on the regular parts, Gre...     

光滑扩展有限元方法及其特点研究

将光滑有限元法S-FEM（Smoothed Finite Element Method）的子域光滑应变技术和边域光滑应变技术同时引入到扩展有限元XFEM（Extended Finite Element Method）中，提出一种新的光滑扩展有限元法S-XFEM（Smoothed Extended Finite Element Method）。在单元选取及扩充结点选取时采用ES-FEM的光滑域划分方式，在数值积分计算刚度矩阵时采用基于三角形子域的CS-FEM积分思路，并给出了高斯点的积分策略。设计了S-XFEM程序架构并利用Matlab语言编制了S-XFEM计算程序。通过几个经典算例研究对比了XFEM和S-XFEM的特点，验证了S-XFEM的精确性和适用性。结果表明，XFEM和S-XFEM均具有很高的计算精确性和收敛性，XFEM计算精度略高于S-XFEM，而S-XFEM在网格独立性上则明显优于XFEM。     

基于一类非局部宏-微观损伤模型的裂纹模拟

结合近场动力学和统一相场理论的基本思想, 最近提出了一类非局部宏-微观损伤模型, 为固体裂纹扩展模拟提供了新途径. 本文在此基础上改进了微观损伤准则, 并给出损伤的$\bar{\lambda}-\ell$语言以刻画固体破坏过程中位移场的不连续程度. 在改进模型中, 首先根据两物质点(即物质点对)之间的变形量, 基于相对临界伸长量的历史最大超越程度, 给出表征物质键性能退化的微细观损伤. 进而, 对影响域内的物质键损伤进行空间局部加权平均, 获得宏观拓扑损伤. 通过引入能量退化函数, 建立基于能量的损伤与宏观拓扑损伤之间的关系, 由此将其嵌入连续损伤力学基本框架, 形成了问题求解的基本方程. 该模型是一类非局部化模型, 可采用有限单元法进行离散求解, 避免了经典局部损伤力学所面临的网格敏感性问题. 文中, 进一步将其应用于具有强非线性回弹特性的裂纹扩展模拟问题. 实例分析表明, 本文方法不仅可以把握裂纹扩展模式, 而且能够定量刻画裂纹扩展过程中的载荷-变形关系. 最后指出了需要进一步研究的问题.      

Characterization and modeling of slow crack growth behaviors of defective high-density polyethylene pipes using stiff-constant K specimen

In this study, slow crack growth (SCG) resistances of defective and normal high density polyethylene (HDPE) pipes were measured using the stiff-constant K (SCK) specimen, where the stress intensity factor (SIF) was maintained at a constant value within a certain crack length range. A significantly reduced SCG resistance was observed in the defective pipe; a detailed procedure for evaluating SCG kinetics using the SCK specimen has been provided herein. The results of a fracture surface analysis indicate that the white window patterns, resulting from poor carbon black dispersion, are the main reason for poor SCG performance. In addition, a crack layer (CL) model was derived for the SCK specimen geometry and was compared with experimental results. It was observed that the crack and process zone growth resistance parameters were significantly lower in the case of the defected pipe than those in the case of the normal pipe.     

Effects of fabric architectures on mechanical and damage behaviors in carbon/epoxy woven composites under multiaxial stress states

The mechanical properties and damage evolutions of carbon/epoxy woven fabric composites with three different fabric architectures, including one plain weave and two twill weave patterns, are experimentally investigated under multiaxial stress states. In particular, the effects of weave patterns are investigated by monotonic and cyclic off-axis tension tests. Both elastic modulus and strength degrade remarkably with increasing off-axis loading angle, while Poisson's ratio is much higher than that measured from on-axis tests and increases with loading strain gradually. Different fabric architectures show limited effects on the modulus and strength under multiaxial stress states, and they are well predicted by transformation equation and Tsai-Wu failure criteria, respectively. However, significantly different failure behaviors are observed in three fabric composites, and microstructure observation shows that fabric architecture affects the stress concentration and the damage development. Smaller crimp ratio and compacted structure postpone the damage development but result in more abrupt failure under multiaxial stress states.     

Interplay between polarization rotation and crack propagation in PMN-PT relaxor single crystals

Investigations on the interconnection between the polarization rotation and crack propagation are performed for [110]-oriented 74Pb(Mg_1/3Nb_2/3)O₃-26PbTiO₃ relaxor ferroelectric single crystal under electric loadings along [001] direction. The crystal is of predominantly monoclinic M_A phase with scatter distributed rhombohedral (R) phase under a moderate poling field of 900 V/mm in [001] direction. With magnitude of 800 V/mm, a through thickness crack is initiated near the electrode by electric cycling. Static electric loadings is then imposed to the single crystal. As the applied static electric field increases, domain switching in the monoclinic M_A phase and phase transition from M_A to R phase occur near the crack. The results indicate that the crack features a conducting one. Whether domain switching or phase transition occurs depends on the intensity of the electric field component that is perpendicular to the applied electric field.     

SOLUTION AND ANALYSIS OF CIRCULAR TUNNEL FOR THE STRAIN-SOFTENING ROCK MASSES CONSIDERING THE AXIAL IN SITU STRESS AND SEEPAGE FORCE

为了研究轴向应力和渗透力共同作用下软化围岩的应力与位移的变化及分布规律. 基于摩尔-库伦屈服准则及应力-应变软化模型并考虑轴向应力和渗透力的共同作用,将整个塑性区分为有限个同心圆环,以弹塑性交界面处的应力、应变为初始值,并采用微小径向应力增量逐步求出各个圆环上的应力应变及塑性区半径,据此重构了考虑渗透力和轴向力共同作用下软化围岩应力应变特性的逐步求解方法. 利用该方法,推导出软化围岩应力应变的解. 计算结果表明：在考虑轴向应力作用下,塑性区半径和隧道围岩位移都随着渗透力的增加而有所减小;当轴向应力为最小主应力时,渗透力的影响更为显著. 这说明渗透力的存在对于隧道围岩的应力应变分布以及塑性半径和围岩的位移有不可忽略的影响.     

Numerical study of contact forces for crack closure analysis

Plasticity induced crack closure (PICC) has been widely studied using numerical models. Different numerical parameters can be considered to quantify the opening level, namely one based on the analysis of contact stresses at minimum load. A modified version of this parameter is proposed here, based on nodal contact forces instead of contact stresses. The predictions were found to be similar to those obtained from the contact status of 2nd node behind crack tip. The PICC_contact parameter was also found to be very consistent and adequate for parametric studies of the influence of different physical parameters. The contributions to the opening stress intensity factor of different points along crack flank were found to strongly decrease with distance to crack tip. The cumulative K_open between the crack tip and a distance of 0.1 mm was found to vary from 30% to 100%, increasing with stress ratio, R. Finally, a K solution was developed for punctual forces applied on crack flank and compared with a literature solution for infinite plates. A good agreement was found for plane strain state but significant differences of about 10% were found for plane stress state.     

Weight function theory and variational formulations for three-dimensional plane elastic cracks advancing

The weight function theory for three-dimensional elastic crack analysis received great attention after the work of Rice (1985, 1989). Several applications have been considered since then, particularly in the context of configurational stability, crack path prediction, stress intensity factor expansions, perturbation approaches. In all cases, a specific hypothesis has been made on the variation of crack shape, in order to formulate the problem in terms of Cauchy principal value. In the present note, such hypothesis is further investigated and consequences discussed. A variational statement given in Salvadori and Fantoni (2013a) is thus rephrased in terms of weight functions. Its discrete formulation shows the potential to accurate approximation of crack front propagation.     

Symmetry breaking in an initially curved pre-stressed micro beam loaded by a distributed electrostatic force

The symmetric and asymmetric buckling of an initially curved micro beam subjected to an axial pre-stressing load and transversal distributed electrostatic force is studied. The analysis is based on a reduced order (RO) model resulting from the Galerkin decomposition with buckling modes of a straight beam used as the base functions. The criteria of symmetric limit point buckling and of non-symmetric bifurcation are derived in terms of the geometric parameters of the beam and the axial load. Two symmetry breaking conditions, defining the relations between the axial load and the geometric parameters of beams for which an asymmetric response bifurcates from the symmetric one, are obtained. The necessary criterion establishes the conditions for the appearance of bifurcation points on the unstable branch of the symmetric response curve; the sufficient criterion assures a realistic asymmetric buckling bifurcating from the stable branches of the symmetric response curve. A comparison between the RO model results and those obtained by direct numerical analysis shows good agreement between the two and indicates that the obtained criteria can be used to predict symmetric and non-symmetric buckling in electrostatically actuated curved pre-stressed micro beams. It is shown that while the symmetry breaking conditions are affected by the nonlinearity of the electrostatic force, its influence is less pronounced than in the case of the symmetric snap-through criterion. The nature of the latter and the relations between it and the symmetry breaking criteria are found to go through a prominent qualitative change as the initial distance between the beam and the electrode, characterizing the electrostatic force, changes.